Effects of thermal motion on electromagnetically induced absorption 
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We describe the effect of thermal motion and buffer-gas collisions on a four-level closed A'' system 
interacting with strong pump(s) and a weak probe. This is the simplest system that experiences 
electromagnetically induced absorption (EIA) due to transfer of coherence via spontaneous emission 
from the excited to ground state. We investigate the influence of Doppler broadening, velocity- 
changing collisions (VCC), and phase-changing collisions (PCC) with a buffer gas on the EIA spec- 
trum of optically active atoms. In addition to exact expressions, we present an approximate solution 
for the probe absorption spectrum, which provides physical insight into the behavior of the EIA 
peak due to VCC, PCC, and wave-vector difference between the pump and probe beams. VCC are 
shown to produce a wide pedestal at the base of the EIA peak, which is scarcely affected by the 
pump-probe angular deviation, whereas the sharp central EIA peak becomes weaker and broader 
due to the residual Doppler-Dicke effect. Using diffusion-like equations for the atomic coherences 
and populations, we construct a spatial-frequency filter for a spatially structured probe beam and 
show that Ramsey narrowing of the EIA peak is obtained for beams of finite width. 

PACS numbers: 42.50.Gy, 32.70.Jz 



I. INTRODUCTION 

The absorption spectrum of a weak probe, interact- 
ing with a pumped nearly-degenerate two-level transi- 
tion, can exhibit either a sharp subnatural dip or peak 
at line center [l[, depending on the degeneracy of the 
levels, the polarizations of the fields, and the absence or 
presence of a weak magnetic field. The phenomenon is 
termed electromagnetically-induced transparency (EIT) 
[3, l3| when there is a dip in the probe spectrum and elec- 
tromagnetically induced absorption (EIA) [J] when there 
is a peak. 

In the case of orthogonal polarizations of the pump and 
probe, both EIT and EIA are related to the ground-level 
Zeeman coherence, which is induced by the simultane- 
ous action of both fields. The simplest model system 
that exhibits EIT is the three-level A system, where the 
two lower states gi^2 are Zeeman sublevels of the ground 
hyperfine level Fg. In a A system, quantum coherence 
can lead to the destructive interference between the two 
possible paths of excitation. As a result, if the pump 
field is tuned to resonance, the narrow dip in the probe 
absorption spectrum at the two-photon resonance can be 
interpreted as EIT caused by a coherent population trap- 
ping in the lower levels. The simplest system that 
exhibits EIA is the four-level A^ system 0, 0] (Fig. [H 
top), consisting of states 511^2 and ei_2 which are Zeeman 
sublevels of the ground [Fg) and excited [F^) hyperfine 
levels, where the gi ^ ei, i = 1,2, transitions interact 
with non-saturating pump(s), and the 32 ^ ei transi- 
tion interacts with a weak probe. The N system gives 
similar results to those obtained for a closed alkali-metal 
Fg ^ Ff. = Fg -\- \ transition interacting with a a+ po- 
larized pump, and a weak tt polarized probe (7|, |8|. It 
has been shown p, 0, Q, that the EIA peak is due to 
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FIG. 1: Top: The A— configuration atom. The light-induced 
transitions are marked by solid (pump field) and dashed 
(probe field) lines, and the wavy arrows are spontaneous decay 
paths. The thick arrow illustrates the spontaneous transfer of 
coherence (TOC). Bottom: probe and pump beam(s), of pos- 
sibly a finite size, propagating through the vapor cell. The 
optical axis is parallel to z, while x and y form the optical 
transverse plane. 



transfer of coherence (TOC) from the excited state to 
the ground state, via spontaneous emission. The excited- 
state coherence only exists in systems where the coherent 
population trapping is incomplete so that there is some 
population in the excited state 0, [l^ . The transfer of 
this coherence to the ground state leads to a peak in the 
contribution of the ground-state two-photon coherence 
to the probe absorption at line center, instead of the dip 
that occurs in its absence (for example, in a A system or 
a non-degenerate A^ system) |ll| . 



In this paper, we investigate the effect of the thermal 
motion of the alkah-metal gas on the EIA spectrum, in 
the presence of a buffer gas. In a previous paper [iJl, we 
discussed the effect of phase-changing colhsions (PCC) 
with the buffer gas on an N system and showed that they 
lead to considerable narrowing of the EIA peak in both 
the presence and absence of Dopplcr broadening. These 
collisions increase the transverse decay rate of the optical 
transitions, resulting in the so-called pressure broadening 
of the optical spectral line, and are thus easily incorpo- 
rated in the Bloch equations. However, in order to de- 
scribe the overall effect of buffer-gas collisions, it is neces- 
sary to include both velocity-changing collisions (VCC) 
as well as PCC [iSl, [ij], which is a much greater chal- 
lenge. Due to the complexity of the problem, we limit 
our discussion to a four-level A^ system, and to buffer- 
gas pressures that are sufficiently low so that collisional 
decoherence of the excited state [1^ can be neglected. 

The Doppler effect occurs in the limit of ballistic 
atomic motion, when the mean free-path between VCC 
is much larger than the radiation wavelength. Due to 
their narrow spectral response, Raman processes such 
as EIT and EIA are much more sensitive to the "resid- 
ual" Doppler effect, arising when there is a difference 
between the wavevectors of the Raman fields. In many 
cases however, the Raman wavelength can become much 
larger than the typical free-path between collisions. For 
example, an angular deviation of a milliradian between 
the two optical beams yields a superposition pattern with 
a wavelength in the order of a millimeter. In this limit, 
the atoms effectively perform a diffusion motion through 
the spatial oscillations of the superposition field, leading 
to the Dicke narrowing of the residual Doppler width. 
While the residual Doppler broadening is linearly propor- 
tional to the Raman wave vector, Dicke narrowing shows a 
quadratic dependence. This behavior was demonstrated 
in EIT with non-collimated pump and probe [1^, [l3 ■ 

Recently, a model describin g th ermal motion and col- 
lisions for EIT was presented [l7H13 |. utilizing the den- 
sity matrix distribution in space and velocity with a 
Boltzmann relaxation formalism. The model describes 
a range of motional phenomena, including Dicke nar- 
rowing, and diffusion in the presence of electromagnetic 
fields and during storage of light. This diffusion model 
was used to describe a spatial frequency filter for a spa- 
tially structured probe [l9| and also Ramsey narrowing 



20|,|21[. Here, we utilize a similar formalism to estimate 



the influence of the atomic thermal motion in a buffer-gas 
environment, including VCC and PCC, on the spectral 
shape of EIA in a four-level N system, with collimated 
or non-collimated light beams. In Sec. [TTl the Dopplcr 
broadening and Dicke narrowing effects are studied for 
plane-wave fields. As the full mathematical treatment is 
lengthy, it appears in Appendix A. However, an approx- 
imate equation which describes the main features of the 
spectra is presented in Sec. [Hi Diffusion-like equations 
for the ground and excited state coherences and popula- 
tions are derived in Appendix B. Two main phenomena 



arc described using this model: (i) a spatial-frequency fil- 
ter for structured probe fields which is presented in Sec. 
mil and (ii) atomic diffusion through a finite-sized beam 
resulting in Ramsey narrowing of the EIA peak, which is 
discussed in Sec. IIVI Finally, conclusions are drawn in 
Sec. El 



II. THE DOPPLER-DICKE LINE SHAPES OF 
EIA 

Consider the near-resonant interaction of a four-state 
atom in an N configuration, depicted in Fig. [TJ The two 
lower states gi and 52 are degenerate and belong to the 
ground level with zero energy, and the excited states ei 
and 62 are degenerate with energy Itljq. The light field 
consists of three beams, each with a carrier frequency ujj 
and wavevector q^, where j = 1, 2 denotes the two strong 
pump beams, and j = p the weak probe, 



E(r,<)= J2 E,(r,t)e- 



-iujjt-^icij r 



C.C. 



(1) 



Here, Ej (r, t) are the slowly varying envelopes in space 
and time. The pumps drive the gi O ei and 32 ■^ 62 
transitions, and the probe is coupled to the 52 <^ £1 
transition. 

Our model will incorporate four relaxation rates: F, 
the spontaneous emission rate from the each of the ex- 
cited states to all the ground states; Fpcc, the pressure 
broadening of the optical transitions resulting from PCC; 
7vcc, the velocity autocorrelation relaxation rate (l/7vcc 
is the time it takes the velocity vector to vary substan- 
tially) [23l, which is proportional to the rate of VCC; 
and 7 is the homogenous decoherence rate within the 
ground and excited state manifolds due, for example, to 
spin-exchange and spin-destruction collisions [2g| . In the 
model, the transition gi -s-)- 62 is forbidden (due to some 
selection rule such as angular momentum). 

To focus the discussion, we assume that all three 
beams arc continuous waves, namely Ej(r, t) = Ej(r). 
We then obtain stationary Rabi frequencies, given by 
Vj — Vj (r) = ^jEj{r)/h, where /ij is the transition 
dipole moment. The complete set of Bloch equations for 
the four-level N system consists of sixteen equations [31 ■ 
In order to simplify the application of the theory to EIA, 
we assume that Vp <g; Vi^2 < P and that the pump tran- 
sitions are well below saturation, so that in the absence 
of the probe, the population concentrates in the §2 state, 
the §2 ■<-> £2 dipole is excited, and the £2 state is empty 
up to second order in the pump field [7|. The equations 
can then be written up to the first order in the probe field 
Vp , which reduces the number of Bloch equations to 
five. 

The complete analytical development is presented in 
Appendix A, and an example of the calculated probe 
absorption spectrum for collincar and degenerate beams 
(qi = q2 = q) is given in Fig. ^a,) (blue line). For the 
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FIG. 2: (a) The probe absorption, calculated from the ex- 
act solution for the density matrix (blue line) and from the 
approximate solution Eq. Q (red line), for collinear plane- 
wave beams (qi = q2 = Qp), Fpcc = 5r, 7vcc = 0.025r , 
A = 0.816, V2 = O.ir, Vi = AV2, Vp = 7 = O.OOir, and 
Ai = A2 — 0. The inset depicts a zoom on the EIA peak 
with its wide pedestal, (b) The EIA spectrum (red line) is 
the sum of three contributions in Eq. (0: the one-photon 
absorption (black dashed line) , a pedestal at the base of the 
peak (brown dotted line), and the sharp peak (green dash- 
dotted line). For the clarity of presentation, the one-photon 
absorption is added to the pedestal and to the sharp peak 
curves. 



numerical calculations, we have considered the D2 line 
of ^^Rb (wavelength 780 nm) at room temperature, with 
a total spontaneous emission rate F = 27r x 6 MHz [23j . 
Other parameters are indicated in the figure caption and 
described in what follows. 

Four complex frequencies control the EIA dynamics, 
each relating to a different coherence in the process: 



Raman (two-photon) detuning. ^4 relates to the three- 
photon transition (whose direct optical-dipole is forbid- 
den), required for the EIA process. Note that the fast 
optical decay rates (F or F) is absent only from ^1. 

In EIA, in contrast to EIT, a strong optical-dipole 
transition ((72 O 62) is excited even in the absence of 
the probe. Its excitation depends on its resonance with 
the pump field, and is thus affected by Doppler broad- 
ening. This leads to velocity-dependent equations even 
in zero-order in the probe field, and introduces the addi- 
tional complex frequency 



^5 = -A2 + q2 • v+i(F + 7vcc), 



(3) 



with the one-photon Doppler shift q2 • v. The overall 
dynamics is thus governed by the five equations (jA9a[) - 
(liSil) . 

We start by calculating the probe absorption spectrum 
for uniform pump and probe fields (plane waves) by solv- 
ing the equations analytically. The spectrum depends on 
18 different integrals over velocity, of the form 



G, = d'v 



-;3„.?" 



•e^ 



S.5U 



^(v), 



(4) 



where -F(v) = (27ru^iJ e~^ Z^""' is the Boltzmann 

velocity distribution, and w^j^ = ki,T/m is the mean ther- 
mal velocity. The determinant ^d, 

U = 666^4 - 6(6 V? + 6V^i') + ^VlV2bAT (6 + 6) , 

(5) 
introduces the power broadening effect (first and second 
terms) , i. e. the dependence of the Raman spectral width 
on the pump powers, and the spontaneous TOG from the 
excited state to the ground state (last term). The last 
term is associated with the TOG due to its dependence 
on the parameter 6, which sets the amount of TOG in 
the original dynamic equations (jAip , and can take either 
the value (no TOG) or 1 Q. The spontaneous decay 
branching ratio is given by A^ ^ y^^^J {^il^g^ + ^^l^g2) 
[31 . The TOG term in Eq. (O depends on the complex 
frequency 



6 = (Ap - Ai) - (qp - qi) • v+i(7 + 7vcc), (2a) 

6=Ap-qp-v-Hi(f -H7vcc), (2b) 

6 = (Ap - A2) - (qp - q2) ■ v+i(F + 7 + 7vcc), (2c) 

6 = (Ap - Ai - A2) - (qp - qi - q2) • v+i(f + 7vcc), 

(2d) 



with the one-photon dctunings Aj = ujj —uje^gj (.7 = 1,2) 
and Ap ~ ujp — uJeig^ , and F = F/2 + Fpcc + 7- The fre- 
quency 6 is related to the probe transition and includes 
the one-photon Doppler shift qp ■ v. 6 and 6 relate to 
the slowly varying ground and excited state coherences 
and include the residual Doppler shift {^p — ^i)■'v and the 



6 + 6 = (2Ap - Ai - A2) - (2qp - qi - q2) • V 

+2i(F/2 + Fpcc + 7 + 7vcc). (6) 

It is important to note that, although each of the individ- 
ual frequencies ^2 and 6 is affected by a Doppler shift 
(either one- or three-photon), the sum 6 + 6 exhibits 
only a residual Doppler shift (assuming nearly collinear 
pumps, qi w q2). Nevertheless the relaxation rate 
(F/2 -f Fpcc + 7 + 7vcc) is the same as that characterizing 
the decay of the optical transitions. As a consequence, 
even when Fpcc is much smaller than the optical Doppler 
width, it plays a significant role in determining the in- 
tensity of the EIA spectrum. This is in contrast to one- 
and two-photon processes (such as EIT), in which Fpcc 



is irrelevant when it is much smaher than the Doppler 
width. It can also be seen that when qi « q2, the var- 
ious residual Doppler shifts are negligible compared to 
the relaxation rates in the determinant S.d, so that £^d is 
only weakly dependent on these shifts. 

Examining the absorption spectrum in Fig. [2][a), we 
observe the narrow absorption peak on top of the broad 
one-photon curve. Moreover, as can be seen in the inset, 
the EIA resonance consists of two independent features: 
a "pedestal" at the base and a sharp absorption peak 
at the center. In order to obtain physical insight into 
these features, we have derived an approximate solution 
for the probe absorption which incorporates the main 
contributions to the EIA, namely the underlying EIT 
mechanism plus the spontaneous TOC. The approximate 
Fourier transform of the nondiagonal density-matrix ele- 
ment for the probe is 
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and rig is the number density of the active atoms. It can 
be shown that Eq. ([7]) is valid provided 7vcc ^ rpcc+r/2. 
For an atom at rest and in the absence of collisions, so 
that 7vcc = 0, Wth — > 0, and Fpcc = 0, Eq. (O is identical 
to the expression obtained by Taichenachev et al. [(| 
(with b = 1), 
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m'/r 



iAr. 



(8) 

The first term in the square brackets in Eqs. ([7]) and ([5]) 
describes the one-photon (background) absorption, and 
the other terms are the EIA peak. 

For a moving atom, the spectrum resulting from Eq. 
([7]) is shown in Fig. [2Ia) (red dashed line) and is com- 
pared with the exact solution; evidently, there is a good 
agreement between the spectra. Despite the small dis- 
crepancy in the intensity of the sharp peak, the approx- 
imate solution preserves the main features in the reso- 
nance. When plotted separately in Fig. [2Ib), the three 
terms in Eq. ([7]) can be identified with the different spec- 
tral features: — G4 (black dashed line) describes the back- 
ground absorption; V2G5 (brown dotted line), which con- 
stitutes the total peak in the absence of VCC, describes 
the wide pedestal; and iG2G37vcc/(l — Gi7vcc) (green 
dashed-dotted line) describes the sharp EIA peak, in- 
duced by VCC. 

Fig. |3] shows the effect of varying the VCC rate, for a 
fixed PCC rate (Fpcc = 5r) and zero pump-probe angular 
deviation. The width of the pedestal feature depends on 
the VCC rate and is given by 7vcc + 7, while the width 
of the narrow peak shows only a very weak dependence 
on 7vcc- Increasing the VCC rate leads to a decrease 
in the overall EIA intensity, but to an increase in the 
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FIG. 3: The EIA peak for different 7vcc rates and Fp 



5r 



at zero pump-probe angular deviation; other parameters as 
in Fig. [3 



ratio between the amplitude of the narrow peak and the 
pedestal baseline. 

We now turn to explore the residual (two-photon and 
four-photon) Doppler and Dickc effects due to wave- 
vector mismatch between the pump fields and the probe, 
introduced in principle either by a frequency detuning 
between the fields, |qp| 7^ |qi,2|, or due to an angular 
deviation between them, q^ j| qi_2- We mainly focus 
on the latter, which may be found in a nearly degen- 
erate level scheme, and we further take the two pump 
fields to be the same, namely qi = q2. Figure [H presents 
the probe absorption spectrum for different values of the 
wave-vector difference, ^q = qp — qi.2, when 7vcc = O.IF 
and Fpcc = F. As can be seen, increasing i5q broad- 
ens the EIA spectrum (see inset). This is analogous to 
the broadening of an EIT transmission peak in a sim- 
ilar configuration [17[. However, the wide collisionally- 
broadcned pedestal remains unaffected by the changes in 
(5q, indicating that it mostly originates from homogenous 
decay processes. Figure [SJ^a) summarizes the full- width 
at half-maximum (FWHM) of the EIA peak for Fpcc = T 
and for various values of 7vcc7 as a function of (5q. Be- 
cause of the difficulty of separating the sharp peak from 
the background in the calculated spectra [231 , the widths 
of the sharp EIA peak were obtained only from the third 
term in Eq. ([7]). In contrast to an EIT peak, which 
does not depend on 7vcc when 5q = (coUinear degener- 
ate beams) [3, the FWHM of the EIA peak at (Sq = 
depends weakly on the VCC rate (although barely no- 
ticeable in the figure). This difference derives from the 
effect of collisions on the pump absorption in the case of 
EIA, as described earlier. 

For (5q 7^ the FWHM of the peak in the Dicke limit 
(high 7vcc) depends on 7vcc and is proportional to the 
residual Doppler-Dicke width, 2i'th^q^/7vcc- In this limit, 
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FIG. 4: Calculated probe absorption spectra with 7v 

and Fpcc = T, for different pump-probe angular deviations. 

Other parameters as in Fig. [2l 



the results are well approximated by the analytic expres- 
sion [ill [dotted lines in Fig. ^a)]: 



FWHM = 2 X — 7,cc 



H 



VthSq 



(9) 



where H (x) 



- 1 



and 



2/ In 2. Increasing 



the pump-probe angular deviation reduces the efficiency 
of the EIA process and thus results in a decrease in the 
probe absorption [Fig. EJb)]. This is of course the oppo- 
site trend to that of EIT (blue stars), where the depth 
of the dip decreases (the absorption increases) with in- 
creasing Sq |19| . 



III. SPATIAL-FREQUENCY FILTER 

We now to turn to discuss the results of our model from 
the viewpoint of a spatial-frequency filter for a structured 
probe beam. When non-uniform beams are considered, 
the different spatial frequencies that comprise the beams 
result in different Doppler and Dicke widths. Conse- 
quently, the various spatial-frequency components expe- 
rience different absorption and refraction in the medium. 
Specifically, the dependence of the absorption on the 
transverse wave-vectors of the probe beam manifests a 
filter for the probe in Fourier space. 

We assume an optical configuration of two coUinear 
uniform pumps (plane waves with Vi and V2 constant) 
and a spatially varying propagating probe, Vp = Vp{r,t). 
Since the medium exhibits a non-local response due the 
atomic motion, the evolution of the probe is more natu- 
rally described in the Fourier space VJ,(k, w) where k and 
Lo are the spatial and temporal frequencies of the envelope 
of the probe. Under these assumptions, the model re- 
sults in a Diffusion-like equations for the populations and 
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FIG. 5; Calculated EIA FWHM (a) and absorption (b) for 
Fpcc ~ F and various values of 7vcc, as a function of the 
pump-probe wave-vector difference 5q. The blue starts are 
the EIT absorption on the Raman resonance conditions, with 
Fpcc = F and 7vcc = 0.025F. 



coherences of the atomic medium, derived in Appendix 
B. To simplify the general dynamics of Eqs. ()B8aP and 
(|B9|) . we take the stationary case [oj = 0, Vp — V^(k)] 
and assume that the carrier wave-vector of the probe 
is the same as that of the pumps, Qp = qi = q2, so 
that (5qi = Sq2 = 0. Taking the Fourier transform [see 
Eq. (jAlip ]. we obtain a set of steady-state equations 
for the spatially-dependent atomic coherences, i?g^g2(k), 



i?eie2(k), and Re^g^(]i.) 



R. 



9192 

(10a) 



i (Ap - Ai) - 7 - ifip \ViY - A'3p 1^21' - Dk" 
= {D^ - bAT)Re,e, + KipV{Vpno, 
[iiAp~A2)~T-j~De]Re,e, 

= -ViV2*{Kip + K3p)Rg,g, - V2*{Kip + Kpump)Vpno, 

(10b) 

i?ei32 = «^'lp {VlRgig2 + VpHo) (lOc) 



where ATip = iGip/ (1 — iGip7vcc) is the one-photon ab- 
sorption spectrum with Gip = J F (v) /£,2(Pv ; K^p = 
iGzpl (1 — iG3p7vcc) is the three-photon absorption spec- 



trum with G. 
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j F{-v) lUd^v; and Kp 



iGpump/ (1 - JGpump7vcc) IS the one-photon (pump) ab- 



sorption spectrum with Gpu 



/F(v)/e5rf-V, as de- 



scribed in Appendix B. Solving Eq. pU]) for i^eiga (k, w), 
substituting the result into the expression for the linear- 
susceptibility [Eq. (JA14P ]. assuming that Vi = 77 V2 
(0 < ?7 < 1), and neglecting all the terms proportional to 
1/r, we obtain 
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FIG. 6: The EIA spatial-frequency filter, given in Eq. (|llb|l . 
as a function of k, with 7vcc ~ 0.025r and Fpcc ~ lOF. Red 
curve is plotted for Ap on resonance, blue-dashed and black- 
dotted curves demonstrate the behavior at nonzero Raman 
detuning. 
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FIG. 7: Calculated probe absorption spectra (blue line) for 
one-dimensional stepwise beam with finite thickness; (a) 2a = 
100 fim and (b) 2a = 10 mm, and fitted Lorentzian (red 
dashed line). All other parameters are the same as in Fig[2] 



Xeiga (k) = -iKuq (1 + L) , 

r;(26^-77)rp 
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7 + 7vcc)] (Pv for 

Ap ^ Fhom = 7 + Tp. In the case where r] = A, Eqs. 
([TT|) is similar to Eq. ([8|) obtained by Taichenachev et 
al. [6|, except for the diffusion term Dk^ , which vanishes 
for an atom at rest. 

The imaginary part of the susceptibihty in Eq. ()11|) 
yields the absorption of the probe for various values of 
k. The first term in the brackets in Eq. piap is the 
linear one-photon absorption, and the second term is the 
fc-dependcnt EIA contribution. Thus, the real part of L 
in Eq. (jllbl) describes an "absorbing" spatial-frequency 
filter, the same way as was done for EIT [2J, [23- Fig. 
ini summarizes several examples of the EIA spatial filter 
behavior as a function of k for Ap = 0, Ap = ±Fhom7 
and Ap = ±2Fhom- At Ap = 0, the curve is a Lorentzian 
and maximum absorption is achieved. When Ap ^ the 
filter becomes more transparent. 



of time in the interaction region and then leave the light 
beams, evolve 'in the dark', and diffuse back inside. Such 
a random periodic motion was described recently by Xiao 
et al. [20, [21| for an EIT system, and was shown to result 
in a cusp-like spectrum. Near its center, the line is much 
narrower than that expected from timc-of-fiight broaden- 
ing and power broadening, and the effect, resulting from 
the contribution of bright-dark-bright atomic trajectories 
of random durations, was named Ramsey narrowing. 

Ramsey-narrowed spectra can be calculated analyti- 
cally from the diffusion equations of the atomic coher- 
ences when the light fields of both the probe and pump 
beams have finite widths [ly • The EIA spectrum result- 
ing from a one-dimensional uniform light-sheet of thick- 
ness 2a in the x— direction is derived analytically in Ap- 
pendix B [Eq. (|B12p ]. In Fig. [3 we show the spectrum 
for two different thicknesses and the fitted Lorentzian 
curves. Near the resonance, the EIA line for the 100 
jim sheet is spectrally sharper than the fitted Lorentzian 
- the characteristic signature of Ramsey narrowing. In 
contrast, the EIA peak calculated for a 10 mm beam is 
well fitted by the Lorentzian. In addition, the EIA con- 
trast deteriorates as the beam becomes narrower, since 
the interaction area decreases and fewer atoms interact 
with the fields. 



IV. RAMSEY NARROWING 

We now consider the A'^ system interacting with 
coUincar probe and pump beams that have finite widths. 
Due to thermal motion, the alkali atoms spend a period 



V. CONCLUSIONS 

In this paper, we extended the theory that describes 
the effect of buffer-gas collisions on three-level A systems 
in an EIT configuration |17H19| to the case of a four- 
level closed N system which is the simplest system that 



experiences EIA due to TOC. Using this formalism, we 
investigated the influence of collisions of optically active 
atoms with a buffer gas on the EIA peak. In addition 
to the exact expressions, we presented an approximate 
solution for the probe absorption spectrum, which pro- 
vides a physical insight into the behavior of the EIA peak 
due to VCC, PCC, and wave-vector difference between 
the pump and probe beams. VCC were shown to pro- 
duce a wide pedestal at the base of the EIA peak; in- 
creasing the pump-probe angular deviation scarcely af- 
fects the pedestal whereas the sharp central EIA peak 
becomes weaker and broader due to the residual Doppler- 
Dicke effect. Using diffusion-like equations for the atomic 
coherences and populations, the spatial-frequency filter 
and the Ramsey-narrowed spectrum were analytically ob- 
tained. 

In extending the description from the A to the N 
schemes, we have considered several elements that are 
likely to be important in other four-level systems. These 
include the diffusion of excited-state coherences and the 
influence of the thermal motion on the optical dipole in 
the absence of the probe. The latter introduces a Doppler 
contribution into the pumping terms and consequently 
affects the power broadening of the narrow resonances. 



Appendix A: Reduced density matrix 



(i) i 

Here, plJ,' is the density-matrix element of the i— th 
atom (one of many identical particles) to the j—th or- 
der in the probe, and apart from Pgjgl ~ 1, pis — 0. 
We also consider the envelopes of the pumps to be con- 
stant in time so that Vi_2 is shorthand for yi_2 (r). The 
wave equation for the probe field is 



1 92 



dt^ 



A-K d 



72 7572 Ep(r,i) = -^—Pe,g,(r,i), (A2) 



dP 



where Peiga (i";*) = Peigs (r, t) e^*"p*e~''ip'* is the con- 
tribution of the ei ■<-> (72 transition to the expectation 
value of the polarization, Pe^g^ is the slowly varying po- 
larization, and V^ is the three-dimensional Laplacian op- 
erator. With Eq. ([T}, and assuming without loss of gen- 
erality that qp = zgp, as shown in Fig. [IJb), Eq. (jA2[) 
can be written in the paraxial approximation as 



d_ 
di 



'dz 



2(7p 



*^ Vi ) Vp (r, t) = i-^Peis. (r, t) , 



(A3) 
where V^ is the transverse Laplacian operator, and g = 

27ra;p l/^eigal J^ is a coupling constant . 

Following [l9| , we introduce a density- matrix distribu- 
tion function in space and velocity. 



Consider the near-resonant interaction of a light field 
consisting of one or two moderately strong pumps and 
a weak probe, as given in Eq. ([1]), with the four- level 
degenerate N system of Fig. [IJa). We use the first-order 
approximation in the probe amplitude, V^^ and assume 
that V2 <T,Vi ^V2tVp < Vi,2- Since the pump transi- 
tions are assumed non-saturated, the atomic population 
in the absence of the probe concentrates in the 172 state, 
and the population in other states can be neglected. The 
g2 ^ 62 dipole, excited in the absence of the probe, is of 
importance and is thus considered. The resulting Bloch 
equations are 






f^eig2 



"eie2 



r9ie2 



(Wp - Wi) = - [i {Ue.g^ - Weisi ) + 7] pi\Yi 



(Ala) 



+ ^Kpl^g:~^V2p';^^: + bArp2^,l, 

(Ujp) = - [lL,,,g, + r/2 + Fpec] pilY^ 

+^ypp^°J9:+'M% (Alb) 

(Wp -U}2)^- [i (Weig^ - We^gJ + T + T] pije^ 

+ ^V,p'^^l + ^\\p^^:~^V^p'^g:, (Ale) 

(cjp - Ui - U2) = - [i (Weiga ^ ^eigi " ^^€292) 

+ F/2 + Fpec] pg^gl - iVrP^g^gl ' (^Id) 



i^i (-^2) = - [r/2 + Fpec - lU^e292\ Pf^el ("^2) 



+ «^2* A 



(0), 



92 92 



(Ale) 



ss' rss 



, (r, V, t) = ^pL, (i) <5 (r - r, (i)) <5 (v - v, (t)) , 

(A4) 
where the time dependence of p\g, (t) is determined by 
Eqs. (jAip . Differentiating Eq. (jA4p with respect to 
time, we arrive at 



d 



d 



Q-^Pss' +--Q-,Ps 






col 



J2jr^Pls' (t) <5 (r - r, (t)) <5 (v - v, (t)) , (A5) 



where the effect of velocity-changing collisions is taken 
in the strong collision limit in the form of a Boltzmann 
relaxation term 123. 



di'- 



col 



p,,, (r,v,t)-4,, (r,t)-F(v) 



(A6) 

with R^^, = R^^, {r,t) = J d^wp^^, {r,v,t) being the 
density-number of atoms per unit volume, near r in 
space, and 



F ^ F(v) = (27rv. 



'th) 



-3/2 vV2i,th 



, wth = (A7) 

m 



is the Boltzmann distribution. 

Before writing the coupled dynamics of the internal 
and motional degrees of freedom, we introduce the slowly 



varying envelopes of the density-matrix elements, pa 
Pas' (r,v,i), as 



Pgig2 ~ Pgia2^ 



-i(Wp-(^l)t i(qp-qi)-r 



rei32 Keig2^ ^ J 



Peie2 — Peie2^ 



(qp-q2)-r 



X — „ p-«(<^p-"i-"2)t„i(qp-qi-q2)- 

^3162 ~ Pgie2^ ^ 

O — O „i'^2i„-iq2T 

rg2&2 — rg2'^2^ ^ ' 



(A8) 



and similarly the slowly varying densities R^^, 
J d^'vpgg, ■ Eqs. (|Aip then become 



d d 

ot or 



PsiS2 7vcc^gig2-^ 
« (^rPeiS2 ^ ^2/59162) + ^^rpe^e^i 

Pci32 ^ ^vcc-t'-eig2-^ 



dt dr 

= i[VpnoF + ViPg,g^] , 

d_ _^_ .. 

9t 9r 



Peie2 '7vcc-'*^eie2^ 
« (^lPsie2 - ^2*PeiS2) + *^pP32e2 7 



dt dr 



1^2 Pgig2J 
d d ' 

ot or 



^0262 7vcc-Kg2e2^ 



(A9a) 
(A9b) 
(A9c) 
(A9d) 
(A9e) 



where ^i (i = 1 ~ 5) are given in Eq. ^. The ex- 
pectation value of the polarization density Peig2 (''lO ii^ 
terms of the number density Re^g^ (r, i) is Peig2 {^^t) = 
^*i32-^ei92 (r, i), and Eq. (|X3)) becomes 



(AlO) 

We now consider the case of stationary plane-wave 
pumps. For this ease, it is convenient to introduce the 
Fourier transform 



W - k ■ V -F ^1] pg^g^ - ij^ccRgig2 {K w) F 

= {V2Pgie2 - V*Peig2) + ibATp^^^^_, (Al2a) 

a; - k • V -f ^2] Peiga - «7vcc-Reig2 (k, w) F 

~{Vpn^F + ViPg,g^), (A12b) 

a; - k • V -F 6] Peie2 - «7vcc-Reie2 (k, w) i^ 

(^2*Pei32 - ^1^3162) - ypPg2e2, (Al2c) 

W - k • V + ^4] /9gje2 - ■ilvccRgie2 (k, w) i^ 

y;P3,,,, (Ai2d) 

W - k • V + ^5] /9g2e2 - «7vcci?32e2 (k, w) i^ 

V^noF, (A12e) 



and Eq. (jAlOp as 



ik;, - i— + i-— \ Vp (k^uj) ^ i-Reig2 (k, w) . (A13) 



The linear susceptibility Xeig2 (k, w) is defined by 



Re^g^ (k, Uj) = XeiS2 (k, w) -Vp (k, w) 



(A14) 



In order to find the probe absorption spectrum, we 
solve Eqs. (jA12p analytically, obtain an expression for 
Pas', and formally integrate it over velocity. This leads 
to an expression for Rsa' in terms of integrals over veloc- 
ity, in the form of Eq. (JH), such as Gi = / d^v^^^^j^^^, 
which can be evaluated numerically. In the general case, 
the resulting expression for Rsa' is very complicated and 
is not reproduced here. In order to explore the under- 
lying physics, we developed an approximate expression 
for the Fourier transform of the density-matrix element 
that refers to the probe transition, namely, Reig2 [see Eq. 
©]■ 

One can verify that in the absence of the pumps 
{Vi = V2 = 0), the resulting one-photon complex spec- 
trum simplifies to the well known result for the strong 
collision regime, K — iG / {1 ~ ijvccG) , where G = 
Jd^vF/{uj-k-v + C2) ^■ 



fir,t) 



+ C)0 +00 

(Pk ■^^^ f did 



^ikr 



2tt 



27r 



e— */(k,c.), (All) 



and write Eqs. (|A9p as 



Appendix B: Diffusion in the presence of fields 

In order to obtain diffusion-like equations for the 
density-matrix elements and the probe fields, we begin 
by integrating Eqs. (jA9al) and (IA9c|) over velocity and 
obtain 



9 .. 
or 



J 



9192 



--i(Aj,-Ai)+7 



i? 



I {V*Re,g, - 1^2i?3ieJ + bATR,,, 



9192 

(Bla) 



_9 
9r 



iS q2 



3e 



- I (Ap - A2 



r + 7 



^6162 — * V ''^1^91 £2 ^^2 ^£192 T" ''^P -'''92 £2 J 5 

(Bib) 



where Js^' = Jss' (r,i) = J (Pwpss' is the envelope of 
the current density. Expanding Pgig2 ^^'^ Peie2 in Eqs. 
(|A9a|) and (^ A9c|) a s p^s' = -Rss'-F+I/TvccpIs/, muhiplying 
Eqs. (|A9a|) and (|A9c|) by v, integrating the resulting 
equations over velocity using 



d d 

(fvjV^—Rss'F ^ S^JVth-g—Rss', (B2) 



defining the current density of the density matrix by 



In order to calculate Re 



^9262! and J 6^32 1 



"-ei92J ^gie2 ^ 

Jgiea, Jff2e2J we assume in Eqs. (|A9bp . (|A9d|) . and 
(jA9e| that the envelopes change slowly enough such that 
\d/dt + V • d/dr\ < |6,4,5| , and get 



'i^2Peig2 ^IvccReigiF + l [VpUoF + VlPg^g^) , (B6a) 
-i^iPaie2 =7vcci?gie2-F'-«^2*P9l92: (B6b) 

-«?5Pff2e2 ^lwccRg2e2F -iV2noF. (B6c) 



Solving Eq. ()B6| formally for Peigai Pgie2-i ^9262 and sub- 
stituting only their leading parts, i.e. pss' — Rss'F, we 
find 



Peig2 = [lvccReig2 " "^^1-^9192) " ^pUo] F/£,2, (B7a) 



^2*^9192] ^/'?4, 



Pg\<^2 — Yl^cc^g\e2 

^9262 = [7vcci?S2e2 + ^2*"-o] 'F'/Cs 



(B7b) 
(B7c) 



Integrating Eqs. (|B7|) over velocity we get 



7vccJss' = / d^vjp[l], 



(B3) 



and retaining the leading terms in l/7vcci we obtain 



R. 



ei92 



iiS'lp [VlRg^g^ + VpUo] 



Rgie2 — *^3p'^^9l92' 
-^^9262 ^~ ^-"-puinp ''^2 ^Oi 



(B8a) 
(B8b) 
(B8c) 



J9192 + D 



dr 



iSqi 



^gii 



bAT . 



(|/l Jeig2 ^2'igig2j •'€162: 

/vcc /vcc 

d 



"eie2 f -^ 



9r 



+ i(5q2 



i?e 



7v 



''l"Jgie2 ''2 "^6132 T" ''pJg2e 



(B4a) 



(B4b) 



where D = Uth/7vcc- Substituting Jgig2i Jeie2 from Eq. 
dm into Eq. ([BT]), we get 



-"-91 £ 



^-*(Ap-Ai) + 7-'D(^+*<^qi 



i iV^Re^g2 ~ ^2^9162) + bArRe,e2 ~^[^+ ^'^^l 



1^1 •^ei52 ^'2*^gie2) •^6162 

■Tvcc Tvcc 

-^^{Ap-A2)+T + -f-D( — + ^S^2 



(B5a) 



i?e 



« (ViRg^ei - ^2*-Reig2 + ^p-R92e2) ~ ^ ( i~ + *'^^2 



a 



(l'lJgie2 ''2 •'ei92 + '4'"^92e2) 



(B5b) 



where K\^ = iGip/ (1 — Gip7vcc) is the one-photon 
absorption spectrum with dp = J F/S,2d^v , K^p = 
iG'i-pl (1 — G3p7vcc) is the three-photon absorption 
spectrum with Gsp = / F/S,a<Pv and i^pump = 
iGpump/ (1 - Gpump7vcc) IS the one-photon (pump) ab- 
sorption spectrum with Gpump = J F/S^^cPv. In the 
case of coUinear pump and probe beams Jq =(5qi_2 = 
qp ^ qi,2 = 3qz, Eqs. (jB5P and (jBSP form a closed set 
when 



a \ zVi,2 (r) , 



dr 
d_ 
dr 
d_ 
dr 



«<5qi^i 



i^q^ 



Tvcc 
Tvcc 

iVpir^ 

7vcc 



can be neglected in Eq. (jB5[) . These terms vanish com- 
pletely in the special case of pump and probe which are 
plane waves {d/dr = 0), and also coUinear and degen- 
erate (Jq = 0) . They can also be neglected whenever 
1^1, 2. pi ^ 7vcc as is the case in many realistic situations. 
However, the term {d/dr + iJqi) • 6^r/7vccJeie2 in Eq. 
(|B5ap cannot be neglected in the case of coUinear pump 
and probe beams since bAT /^^^c does not go to zero. 

Substituting Eq. (|B4b| into Eq. (|B5a|) . and Eq. ((B8)) 
into Eq. (jBSp . we find: 
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^ - * (Ap - Ai) + 7 + i^lp |Vl|' + A-3P |V^2|' } Rg,g, 



= D [— +iSqij Rg^g^+D i — +iSq2J Re^e^ 



d 

— -i(Ap-A2)+r + 7j>i?eie2 



(B9a) 



£> ( ^ + i5ci2 ) Re^e2 + VlV^* (A'lp + Agp) Rg,g., 



V2 {Kip + Apurap) VpUQ. 



(B9b) 



These are the final difli'usioii-fike coupled equations for 
the ground- and excited-state coherences. 

In order to investigate the Ramsey narrowing of the 
EIA peak, we consider finite probe and pump beams 
and restrict the discussion to coUinear EIA. We assume 
that the fields are stationary and overlap in their cross 
sections with negligible variation along the z— direction, 
Vp{v,t) = Vpw{v^), Vi{v) = Viw{ri_), ^2 (r, i) = 
V2w(rj^), where w (r_L) is the transverse profile of the 
fields. We further take 5q = Q and Ai = A2 = for 
brevity. In the diffusion regime, we rewrite Eqs. 
and (Ell) as 



The solution for Rg^g2, symmetric in x and decaying as 

a; I — ?► 00, is given by 



Rg^g._ {\x\ < a) = 

C2 cosh {kix) + Ci cosh (fc2x) 



Mr/32 + Dall^i 

{Daia2f + bATl33' 

(Blla) 



Rete2 (kl < a) 



Ci (fc| - af) D7vcc 
bAT {Dal + 7vcc) 



C2 {kf ~ al) Dj,,, 
bAT {Dal + 7vcc) 

^9132 (kl > a) = 
C3bAT{Dal+j^ 



■ cosh(fcia;) -t- 



cosh(fc2x) 

/3i/33 - P2C^lD 

PsbAT - {Daia2f ' 
(Bllb) 



{al~al)Dlv. 
Re.e, {\x\ >a) = C3e-"^(l^l-") 



^^a2{\x\-a) _^ (j^^^a,i\x\-a) ^ 



(Bllc) 
(Blld) 



where aj = {-lAp + -f + Kip \Vi\' + Agp \V2n/D, al = 
al + r/D, al = {-iAp + 7) /D, and pi = V^VpKipUo, 
p2 = ViVi{Kip + K3p),P3 = V^Vp{Kip + Kp^^p)na. The 
complex diffusion wave-numbers are obtained from 



iAp + -1+[Kip\Vi 
D 



K- 



3p 



\v2n w{v^f 



Rgig2 — 



bAT 1 



R 



7v 



-vi 



Re-,e2 - KlpV^VpUQW (r_L) , 



eiga = iKip {ViRg^g2 + VpUo) w {r±_) 



UK, 



r + j-Dvl)Re 



= Vi{Kip + K3p)Rg,g2V2*w{r^Y 
+ Vp {Kip + Kpump) noV2W {r±f 

^gie2 

Rgie2 = -iKpumpV.^now (rj_) . 



-iK3pV2Rg,g2W{r±). 



(BlOa) 
(BlOb) 



(BlOc) 
(BlOd) 
(BlOe) 



We further consider a probe and pump beams with a 
uniform intensity and phase within a sheet of thickness 
2a in the a;— direction (one-dimensional stepwise beams): 



w{x,y) 



1 for |.t| < a 
for |a;| > a 



2Dj^cckl2 = -D7vcca+ + PabAT 
T [(^7vcc)'«i + (isbAT (4 + 2D-f,,,al + P^bAT)] '/' , 

with a±. = al ± al- The coefficients Ci (z = 1 — 4) 
are obtained from the continuity conditions of Rss' and 
{d/dx) i?ss'at |x| = a. From Eq. (jBlObp one finds 



Reig2 (kl < a) = ''■K2 



V\ {C2 cosh(fcia;) 



-f Ci cosh {k2x) 



bAT 132 + Da'^ 



2P1 



{Daia2) + bAVPa 



VpUo 



(B12) 



and the energy absorption at frequency tUp is finally cal- 
culated from P (A) = {hujp/a)lm J_ dx-Reiga (^) ■ Two 
examples for the resulting spectrum are given in Fig. [T] 
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The fact that the 'inner' decoherence rate 7 is shared 
by both the ground and the excited manifolds, does not 
imply that their total decoherence rate is the same; the 
coherence between the two excited states decays via the 
e ^ g relaxation channels and therefore decays much 
faster than the ground-state coherence. Incorporating dif- 
ferent values of 7 for the ground and excited states does 
not lead to substantial changes in the collision-induced 
phenomena explored here. 

When 7vcc — > 0, the third term in Eq. (0 gradually 
vanishes, and the second term {V2G5) is responsible for 
the EIA peak, as indicated by the brown-dotted line in 
Fig.[2Ib). Its width is limited by homogenous broadening 
mechanisms and determined by 7 and Fpcc. 



